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CONSTRAINTS ON FLAVOR-DEPENDENT LONG RANGE FORCES 
FROM NEUTRINO EXPERIMENTS 
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We study the impact of flavor-dependent long range leptonic forces mediated by the L e — L M or 
L e — L T gauge bosons on the solar neutrino oscillations, when the interaction range Rlr is much 
larger than the Earth-Sun distance. The solar and atmospheric neutrino mass scales do not get 
trivially decoupled in this situation even if 613 is vanishingly small. In addition, for a > 10 — 52 and 
normal hierarchy, resonant enhancement of 613 may give rise to strong energy dependent effects 
on the u e survival probability. A complete three generation analysis of the solar neutrino and 
KamLAND data gives the 3cr limits a efl < 3.4 X 10~ 53 and « eT < 2.5 X 10 -53 when Rlr is much 
smaller than our distance from the galactic center. With larger Rlr, the collective LR potential 
due to all the electrons in the galaxy becomes significant and the constraints on a become stronger 
by upto two orders of magnitude. 
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1. Introduction 

Flavor-dependent long range (LR) leptonic 
forces, like those mediated by the L e — L M or 
L e — L T gauge bosons, 1 constitute a minimal 
extension of the standard model that pre- 
serves its renormalizability. The extra U(l) 
gauge boson, albeit nearly massless, would 
escape direct detection if it couples to the 
matter very weakly. 2 Since these long range 
forces violate the equivalence principle, they 
are strongly constrained by the lunar ranging 
and Eotvos type gravity experiments. 3 ' 4 For 
a range of the Earth-Sun distance or more, 
these experiments imply the 2a bounds a < 
3.4 x 1CP 49 , where a denotes the strength of 
the long range potential. 

The coupling of the solar electrons to the 
L e — Lp gauge bosons generates a long range 
potential for neutrinos, whose value at 
the Earth is 1.3 x 10" n eV(a e/3 /10" 50 ). The 
typical value of Am 2 / E for atmospheric neu- 
trinos is ~ 1CP 12 eV, so even with the strong 
constraints above, LR forces affect atmo- 



spheric neutrino oscillations. This allows one 
to put stronger constraints 2 on the couplings, 
52 and aw < 6.4 x 1CT 52 . 



a eM < 5.5 x 10 



2. LR potential due to the Sun, 
the Earth and the galaxy 

The behavior of from the solar core all 
the way to the Earth (r w 215 r ) and be- 
yond is shown in Fig. 1. It is seen that V% 
dominates over the MSW potential Vcc in- 
side the Sun for a > 10~ 53 . Moreover, it 
does not abruptly go to zero outside the Sun 
like Vcc, but decreases inversely with r. Its 
value at the Earth surface is an order of mag- 
nitude larger than the potential there due to 
the electrons in the Earth, so one can neglect 
the latter. 

When Rlr > RgaX, °ur distance from 
the galactic center, the collective potential 
due to all the electrons in the galaxy may 
become significant. We denote the galactic 
contribution to the potential V e p as 



V e f = b a ef3 Nl s jR° m 



(1) 
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small parameters x = A21/A32 ~ 0.03 and 
sin #13 < 0.2, except in a narrow region 
around y e ^ = V r eM /A 32 « 2/3. The exact 
numerical values for mixing angles and mf 
for different values of a for normal as well as 
inverted hierarchy are shown in Fig. 2. 

°« = 10 " 53 m 2 > m 2 l <* =io- 52 



log 10 (r/r ) 

Fig. 1. Comparison of the MSW potential Vcc and 
the LR potential V e ® due to the solar electrons. The 
(Am 2 /2E) values corresponding to E = 10 MeV are 
also shown. 

where AT e ° gal = 10 12 7V© and R° al = 10 kpc. 
The parameter b takes care of our ignorance 
about the distribution of the baryonic mass 
in our galaxy. With Rlr > Rgai, we ex- 
pect 0.05 < b < 1. The value of b may be 
smaller if Rlr is smaller, b = is equivalent 
to Rlr <C -Rgai- The screening effects 7 are 
negligible over the scale of i? ga i. 

3. Masses, Mixings and 

Resonances of Solar Neutrinos 

The LR potential gives unequal contribu- 
tions to all three flavors simultaneously un- 
like in case of the charged current potential. 
As a result, the inclusion of three generations 
in the solar analysis becomes necessary. 

The appropriate Hamiltonian in the L e — 
case describing the neutrino propagation 
can be written in the flavor basis as 

Hf = R Hq R t + V , (2) 

with V = Diag(F cc +F eM , ~V eil , 0). Here R = 
#23(023)-Ri3(0i3)-Ri2(0i2), we have assumed 
that no CP violation enters the picture. One 
can take H Q = Diag(0, A 2 i, A 32 ) with Ay = 
Amfj / (2E). The antineutrino propagation 
is obtained by the replacement V — > — V. 

The Hamiltonian (2) can be analytically 
diagonalized 5 by keeping terms linear in the 




Fig. 2. The angles and m 2 values in matter for solar 
neutrinos for E = 10 MeV, in the case Rlr <C fi ga i ■ 
The m 2 values are correct up to an additive constant, 
so that only their relative values have a physical sig- 
nificance. 



3.1. For a < 10~ 52 

Both 823 as well as 6*13 get only small cor- 
rections. In the limit #i 3m — > 0, the third 
mass eigenstate decouples and the scenario 
reduces to 2v mixing. However, the effective 
matter potential is 5 

V12 w Vcc + "Ml + cos 2 9 23m ) , (3) 

and not V cc + 2V efi as would have been taken 
in a naive 2-gcncration analysis. 

The MSW resonance takes place with 
the modified potential Vi2- For a > 10~ 53 , 
the V12 contribution dominates over V cc and 
the resonance is shifted outside the Sun 
where its behavior is solely determined by 
the LR potential. 
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If ft(-E) is the jump probability at the 
v\ m -vim resonance, the net survival proba- 
bility of v e is 

Pee{E) — (1 — Pi) cf 3P C 12 p cf 3E C 12E 

I 75 2 2 2 2 
"T" C 13 p S 12 p C 13 p C 12 p 

+ (1 — ft) C 13P S 12P C 13P S 12B 
I f) 2 2 2 2 

t rx c 13P c 12 p c 13 p s 12B 

+ S 13P S 13P ■ (4) 

Here ftyp and fty^ are the values of 9ij m 
at the neutrino production point and at the 
Earth respectively. 

3.2. For £*>10- 52 

In this range of a, the V\ m -v 2m resonance 
is always outside the Sun and adiabatic. In 
addition the angle 6>i 3m gets resonantly en- 
hanced when y efl w 2/3 (normal hierar- 
chy). The 0i3 m enhancement corresponds to 
the v 2rn -v 3m level crossing, with an effective 
potential 5 

V 23 = V cc + V efi (l + sin 2 9 23m ) ■ (5) 

The net survival probability of v e is 

Pee{E) — Ci 3 p (?\2P C 13£; C 2 2E 

+ (1 — Eh) c 2 3P s\ 2P c 2 3E s 2 2B 

+ ft S 13P C 13P S 12P 

+ (1 - ft) S 13P S 13P 

+ ft C 13P S 12P S 13P > (6) 

where Pjj{E) is the probability that v 2m and 
vzm convert into each other after traversing 
through this resonance. 

In general Ph ~ at high values of #13. 
For 9i3 < 0.1°, the value of Ph becomes sig- 
nificant. In the range where 0.1 < Ph < 0.9 
(the semi-adiabatic range) , Ph is also highly 
energy dependent. 

4. Constraints from solar 
neutrinos and KamLAND 

We perform a global fit to the data from so- 
lar experiments and KamLAND, using the 
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Fig. 3. Ax 2 values and limits for the L e — as 
well as L e — L T symmetry, with #13 = 0° .The case 
Rlr ^ ^gal i s represented by b = and higher b val- 
ues correspond to larger contributions from galactic 
electrons (see Sec. 5). 

X 2 minimization technique with covariance 
approach for the errors. 8 

The best fit values for the solar param- 
eters are always observed to lie in the LMA 
range 6 with vanishing a ep giving the best fit. 
For a < 10 -52 , the value of x 2 is minimum 
for 9\ 3 = 0°. When a > 10~ 52 , a strong en- 
ergy dependence in the survival probability 
is introduced for 9 < 0.1° through Pjj(E), so 
that the x 2 values for extremely low 9\ 3 val- 
ues become large. However, the region a > 
10 -52 turns out to be excluded to more than 
3<7. Therefore we quote the most conserva- 
tive upper bounds on a, by taking 9\ 3 = 0°. 
These limits are shown in Fig. 3: the 3cr limit 
corresponding to the one-parameter fit in the 
L e — case is a e ^ < 3.4 x 10~ 53 . The corre- 
sponding L e — L T limit is a eT < 2.5 x 10~ 53 . 
The bounds are independent of the neutrino 
mass hierarchy. 

5. LR potential from the galaxy 

When Rlr ~> -R ga i, the net potential V efl = 
V® + V e f l is shown in Fig. 4. For V e f l » 
Am%/(2E), there is no MSW resonance that 
is essential for a good fit to the solar neutrino 
data. Even for lower values of b and a, if V^ 1 
dominates over Vcc at the MSW resonance, 
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the resonance tends to be adiabatic even for 
low energies, so the radiochemical data dis- 
favors the solution. 




Iog 10 (r/r ) 



Fig. 4. The net potential V eM 
ious values of b and a. 



The A\ 2 values as a function of a for 
different b values are shown in Fig. 3. The 

case are 



3<7 constraints for the L e 
a efl < 2.9 x 1(T 54 = 0.1) and a ef2 < 
2.6 x 1(T 55 (6 = 1). For L e - L T , the con- 
straints are a eT < 2.3 x 10~ 54 (b = 0.1) and 
a eT < 2.1 x 10~ 55 (6 = 1). We expect b < 1 
even with generous estimates. The most con- 
servative constraints are clearly with 6 = 0, 
as calculated in Sec. 4. 



6. Concluding remarks 

The long range forces mediated by L e — L 
vector gauge bosons can give rise to non- 
trivial 3-v mixing effects inside the Sun, and 
affect the MSW resonance picture. The an- 
gle 6\3 may even get resonantly amplified 
if a > 10~ 52 for normal hierarchy. These 
effects allow one to put constraints on the 
coupling of the LR forces from the solar and 
KamLAND data. The bounds obtained are 
orders of magnitude better than those ob- 
tained earlier from the gravity experiments 
and atmospheric neutrino data. If the range 
of the forces is larger than i? ga i, the bounds 
become still stronger by upto two orders of 
magnitude. 



A recent paper 9 also has given compara- 
ble bounds on the LR couplings. However, 
they assume one mass scale dominance, ne- 
glecting the effect of the third neutrino alto- 
gether. Moreover, they have not taken the 
galactic contribution into account even when 
the range of the force is more than i? ga i- 
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